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■ We study effects of the pfiysical realization of quantum computers on their logical operation. 

' Through simulation of physical models of quantum computer hardware, we analyze the difflculties 

$_( , that are encountered in programming physical realizations of quantum computers. Examples of 

logically identical implementations of the controUed-NOT operation and Grover's database search 
I algorithm are used to demonstrate that the results of a quantum computation are unstable with 

respect to the physical realization of the quantum computer We discuss the origin of these insta- 
bilities and discuss possibilities to overcome this, for practical purposes, fundamental limitation of 
quantum computers. 

PACS numbers: 03.67.Lx, 05.30.-d, 89.80.+h, 02.70Lq 



00 



Oh 



> 
in 
00 
o 

Tj- ' I. INTRODUCTION 

o ■ 

Recent theoretical work has shown that a Quantum Computer (QC) has the potential of solving certain com- 
putationally hard problems such as factoring integers 0] and searching databases much faster than a conventional 
computer . In most theortical work the operation of a QC is described in terms of highly idealized transformations 
on the qubits The impact of the physical implementation of a QC on its computational efficiency is largely 

unexplored. 

The logical operation of conventional digital circuits does not depend on their hardware implementation (e.g. 
' semiconductors, relays, vacuum tubes, etc.). Dissipative processes suppress the effects of the internal, non-ideal 
(chaotic) dynamics and drive the circuits into regions of stable operation. Conventional digital computers, built 
IV from these digital circuits, are in one particular state at a time and are able to perform logical operations that do 
. ^ [ not depend on their hardware implementation. From the point of view of programming the computer this is very 
■ important. Implementations of algorithms designed to run on a conventional computer will give results that do not 
?H ] depend on the hardware used to build the computer. 
. . . 1 A QC differs from a conventional digital computer in many respects. A QC exploits the fact that a quantum system 
can be in a superposition of states. Interference of these states allows exponentially many computations to be done in 
parallel The presence of the superposition of states is a direct manifestation of the internal quantum dynamics 

of the elementary units of the QC, the qubits. In an ideal QC the qubits are assumed to be ideal two-state quantum 
systems. Therefore, the operation of an ideal QC does not depend on the intrinsic dynamics of its qubits. 

A physically realizable QC is a many-body system in which the quantum dynamics of the qubits is essential to its 
operation. Manipulation of one qubit may cause unwanted motion of other qubits. It is difficult to suppress these 
effects by dissipation because in contrast to the case of conventional digital circuits, dissipation processes have a 
devastating effect on the coherent quantum dynamical motion of the qubits. Therefore a quantum algorithm (QA) 
may yield quantum computation results that depend on the specific physical realization of the QC. Although QA's 
can be designed independent of the QC hardware, the implementation of a QA on a physical realization of a QC (i.e. 
the programming of the QC) very much depends on the hardware of which the QC is built from. We refer to the 
difficulty of programming QC's as the Quantum Programming Problem (QPP). 

Due to the QPP it may be very difficult to develop a non-trivial quantum program for a physical realization of a 
QC. Moreover, there is no guarantee that an implementation of a QA, that works well on one QC will also perform 
well on other physical realization's of QC's. As mentioned above, there are several factors that contribute to the QPP: 
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1) Differences between the theoretically perfect and physically realizable one- and two-qubit operations; qubits 
cannot be kept still during the time that other qubits are being addressed; precision needed to implement 
operations on the qubits. 

2) The effect of coupling of the qubits to other degrees of freedom (dissipation, decoherence) . 

In most theoretical work on QC's and QA's one considers theoretically ideal (but physically unrealizable) 

QC's. Then the QPP is not an issue. The QPP is also fundamentally different from the error propagation previously 
studied in QA's implemented on ideal QC's |p|-pl| since the QPP is due to the specific realization of the QC and leads 
to systematic instead of random errors. 

How does a QPP reveal itself? Consider two logically independent operations (Oi and O2) of the machine. On a 
conventional computer or ideal QC, the order in which we execute these two mutually independent instructions does 
not matter: O1O2 — O2O1. However, on a physically realizable QC sometimes the order does matter, even if there 
are no logical dependencies in these two program steps. In some cases, due to practical problems in manipulating 
individual qubits O1O2 7^ O2O1 and the QC may give wrong answers. Note the qualifier sometimes. There seems 
to be no general rule to decide beforehand which operation and at what stage of the QA the QPP leads to incorrect 
results. At present the only way to find out seems to be to actually carry out the calculations and check the results. 

In this paper we study the relation between the physical realization of QC's and their logical operation. We 
investigate various aspects of the QPP by simulating QC hardware. In this work we only consider effects of the 
intrinsic quantum dynamics of the QC (item 1, see above). The study of the effect of the coupling of the qubits to 
other degrees of freedom (item 2, see above) is left for future research. We demonstrate that the programming of a 
physical, non-ideal implementation of a QC is difficult, even if the QC consists of only two qubits. Herman et al. p^ ] 
investigated the influence of the Ising spin interaction on the quantum dynamics of NMR systems. Although they 
did not address the QPP, the work is similar in spirit to and the present paper as it explores the consequences 
of the difference between the ideal QC's and physical realizations of QC's. For the physical systems and time-scales 
considered in this paper, the effects of the interactions between the spins are negligible. As far as we know no 
experimental data has been published that specifically addresses this, for potential applications, very important and 
intrinsic problem of programming QC's. However, with the QC hardware currently available a test of correct quantum 
computation on a physical realization of a QC is definitly within reach. In this paper we propose two simple QA's 
that may be used to study the QPP in physical realizations of QC's. We also discuss methods to enlarge the region(s) 
of reliable operation. 

The paper is organized as follows: In Section || we describe a physical model of a QC. Our choice is largely inspired 
by NMR-QC experiments [p"^-p^, mainly because other candidate technologies for building QC's [p0[-p4| are not yet 
developed to the point that they can execute computationally non-trivial QA's. As the basic example of a QA we 



take the ControUed-NOT (CNOT) gate In Section III we discuss the implementation of the CNOT gate on 

an ideal two-qubit QC and describe simple, non-trivial QA^ that consist of repetitions of CNOT operations. As an 
illustration of the general nature of the QPP, we also consider a more complicated example, namely Grover's QA to 
search for an item in a database In Section IV we derive the conditions for which the physical two-qubit QC 

we 



will exhibit ideal QC behavior and discuss the generalization of these ideas to n-qubit QC's. Also in Section 
describe the implementation of the QA's discussed in Section IH on a physical realization of a QC. In Section^ we 



demonstrate and analyze the QPP by simulating the time-evolution of (= execute QA's on) the physical model of the 



QC by solving the time-dependent Schrodinger equation. In Section VI we summarize our findings 



II. PHYSICAL MODEL OF A QUANTUM COMPUTER 

Generic QC hardware can be modeled in terms of quantum spins (qubits) that evolve in time according to the 
time-dependent Schrodinger equation (TDSE) 

z^^\m) = H{tmt)), (1) 

in units such that h = 1. For present purposes it is sufficient to consider two-qubit QC's only. The state 

mt)) = «a, i; ^)l i, i> + «(T, i; t)\ T, i> + T; *)! i, T> + «(T. T; t)\ T, T), (2) 

describes the state of the QC at time t. The complex coefficients a(J,, l;t), . . . , a(t, T; t) completely specify the state 
of the quantum system. In the absence of interactions with other degrees of freedom this spin-1/2 system can be 
modeled by the time-dependent Hamiltonian 
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H{t) = -JSIS^ - hlSl - h^S'^ - hlSl - h^S^ - h\Sl - hlSl 

-{hlSl + h^S^) sm{Lut + (j)^) - {h\S\ + hlSl) sm{Lut + (j)y), (3) 

where S", a = x,y, z denotes the a-th component of the spin-1/2 operator representing the j-th qubit, J determines 

the strength of the interaction between the two qubits, and represent the strength of the apphed static (magnetic) 

and apphed sinusoidal field (SF) acting on the j-th spin respectively. For a physical system, = l^i ^^'^ ^2 — 7^1; 
for a = x^y^z where 7 is a constant. The frequency and the phase of the SF are denoted by lo and 0ct. As the Ising 
model, i.e. the first term of (^), is known to be a universal QC [ ^6||37| ], model (||) is sufficiently general to serve as a 
physical model for a generic QC at zero temperature. In terms of spin matrices, the operator Qj measuring the state 
of qubit J is given by 



= 2 - ^1- (4) 

For numerical purposes it is necessary to fix as many model parameters as possible. We have chosen to simulate 
the two nuclear spins of the and "'^'^C atoms in a carbon-13 labeled chloroform, a molecule that has been used in 
NMR-QC experiments 1110. In these experiments /if/27r « 500MHz, /i|/27r « 125MHz, and J/27r w -215Hz @. 
In the following we will use model parameters rescaled with respect to /ij/27r, i.e we put 

J = -0.43xlG-^ hl = l, h^ = 0.25. (5) 

Note that there is a difference of many orders of magnitude between the interaction J and the fields h^. If the 
duration of the SF-pulses is much shorter than 2tt/\J\, the effects of J on the time evolution during these pulses are 
very small. Our numerical experiments (see below) are all performed under this condition. We will only consider 
QC's at zero temperature without coupling to the environment. In this sense we simulate highly idealized NMR 
experiments on a closed quantum system at zero temperature. This allows us to study a concrete physical realization 
of a QC and at the same time focus on the intrinsic quantum dynamics of the QC. 

A QA for QC model (||) consists of a sequence of elementary operations (EO) that change the state j^") of the 
quantum processor according to the TDSE, i.e. by (a product of) unitary tranformations. Each EO transforms the 
input state |^'(t)) into the output state \'^{t + t)) where r denotes the execution time of the EO. The action of an 
EO on the state l^t) of the quantum processor is defined by specifying how long it acts (i.e. the time interval r during 
which it is active), and the values of J and all h's. During the execution of an EO the values of J and all h's are kept 
fixed. 

The time evolution of quantum model (3) is obtained by solving TDSE (^ for model (||). The simulations have been 
carried out with a software tool called Quantum Computer Emulator (QCE) The QCE software simulates physical 
models of QC hardware by a Suzuki product-formula p9| , ^0[ |, i.e. in terms of elementary unitary operations [p]-^. 
For all practical purposes, the numerical results obtained by this technique are exact. A detailed description of the 
QCE software tool can be found elsewhere Q]. 



III. IDEAL QUANTUM COMPUTER 



A. Single-qubit operations 



One qubit or one spin-1/2 system is a two-state quantum system. The two basis states spanning the Hilbert space 
are denoted by | t) = |0) and | j) = |1). Rotations of spin j about 7r/2 around the x and y-axis are basic QC 
operations. We will denote them by Xj and Yj respectively. In matrix notation, they are given by 



V2 I * 1 



(6) 



and 



1 
72 



1 1 

-1 1 



(7) 
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Clearly operations such as and (Q) can be implemented in terms of the time evolution of model ^ by a proper 
choice of the model parameters. Writing \a) = ao|00) +ai|10) + a2|01) +a3|ll) with \b1b2) = |6i)|62) and 6^ = 0, 1 we 
have 




1 i 0\ / ao 

i 1 

V2 1 1 z 

i 1 / V 



(8) 



For example, -^i|ll) = (|11) + i|01))/\/2. Using the same labeling of the basis states as in (ph we have 




(9) 



e.g. ^2111) = (|10) + |11))/V2- The matrix expressions for the inverse of the rotations Xi and denoted by Xi 
and X2 respectively, are obtained by taking the hermitian conjugates of the matrices in (^) and (||). For example, 
F2|ll) = (|ll)-|10))/^/2. 



B. Two-qubit operations: CNOT gate 

Computation requires some form of communication between the qubits. A basic two-qubit operation is provided 
by the CNOT gate. The CNOT gate flips the second spin if the first spin is in the down state, i.e. the first qubit acts 
as a control qubit for the second one, see Table |. The procedure that we use to construct the CNOT gate may seem 
a little ad hoc and indeed to considerable extent it is. There is no unique method to construct QC gates. 

On an ideal QC the CNOT gate can be implemented by a combination of single-qubit operations and a two-qubit 
phase shift operation P defined by the matrix 



P = 



( e^'>"' <d\ 

e'"^! 

6'"^=^ 

\ e'-^^ ) 



Assume that the QC is in a state 

I*) = aolOO) + ai|10) + a2|01) + aalll). 
First we apply to the rotation 12, as defined in (^). This gives 

^2!*) = ^[(ao + 02)100) + (ai + a3)|10) + (02 - ao)|01) + (03 - ai)|ll)]. 

Next we apply to I2I*) the phase shift P 

PY2\^) = ^ [e^^VolOO) + e^^^cillO) -t- e''^^C2|01) + e''^^C3|ll)] , 
v2 



(10) 



(11) 



(12) 



(13) 



where cq = + a2, ci — ai + a^, C2 — 0,2 — oq and C3 =03 — 01. Finally we apply the inverse of the rotation Y2 
F2PF2I*) = ^[(e'*"co - e'*^C2)|00) + (e^^^ci - e'*^C3)|10) + (e'^Vo + e*"*^C2)|01) + (e^-^^ci -f e*^^C3)|ll)] . (14) 

We now determine the angles (pi such that the sequence (14) performs the CNOT operation. Since the CNOT gate 
will not change oq and 02 (see Table ^ we can choose 0o = 4'2- This gives 

Y2PY2\'^) = e'^" [oolOO) + 02101) + e'^^{ai cosa + ias sina)|10) + e'^ia^ cosa + iai sina)|ll)] , (15) 
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where P = a + (p^ — (j)Q and a = {(j)2 — 03)/2. The global phase factor e^^° is physically irrelevant. 

The simplest way to implement the phase shift P is to use the time evolution, i.e. P = e~"^' , of the Ising model 



Hj — — JS1S2 — hSi — hS2, 



(16) 

where the external fields acting on both spins are the same. From (|l6|) it follows immediately that (po = t{J/A + h), 
(pi = 4'2 = — tJ/4 and 03 = t{J/4: — h). Taking into account our choice 0o = 4>2, (f5[) becomes 



Y2PY2YI!) = e^"/2[ao|00) +a2|01) +e-*"(aicosa + m3sina)|10) +e-*"(a3Cosa + misina)|ll) 
Using the same labeling of states as in (ph we have 



Y2PY2 = e'"/2 



f I 

e~'"cosa ie^*"sinQ; 

10 

\0 ie~*"sina e'^cosa 



(17) 



(18) 



Comparing the truth table of the CNOT gate (see Table with the matrix in ( [iq), it is clear that putting a — tt/2 
will do the job (upto an irrelevant global phase factor). In terms of Hamiltonian (|16[), — tJ = tt and h = — J/2. The 
sequence 



CNOT = Y2IY2 



10 
1 
10 
10 



(19) 



performs the CNOT operation on qubit 2 with qubit 1 acting as the control variable. Here we introduced the symbol 
/ to represent the time evolution e^"^' with r = — tt/ J. 



C. Quantum Algorithms 

Any QA can be written as a sequence of the one- and two-qubit operations discussed above. As a simple example 
of a QA we take {CNOTf. On an ideal QC, CNOT^ is the identity operation and hence {CNOTf = CNOT but 
on a physical QC this is not always the case, see below. To illustrate the dependence of the quantum computation 
on the physical implementation and on the choice of the input state we consider two QA's, QAi and QA2, defined by 

Q^i|6i&2) = {CNOTf \b1h2), (20) 
QA2\singlet) = Yi{CNOTf\singlet), (21) 

where \singlet) = (|01) - |10))/y2. We have 

(|01) - |11))/V2 = {CNOTf \3inglet), (22) 

and hence {singlet\{CNOTfQi{CNOTf\singlet) = 1/2. We can obtain a clear-cut answer in terms of expectation 
values of the qubits by applying a 7r/2 rotation of spin 1 

\ll) =Yi{C NOT f\singlet). (23) 

Therefore in (pi]), the CNOT operations are followed by a 7r/2 rotation of spin 1. 

As a more complicated example of a QA, we consider Grover's database search algorithm to find the needle in a 
haystack. On a conventional computer, finding an item out of of N elements requires 0{N) queries Grover has 
shown that a QC can find the item using only 0{\^) attempts |Q. The reduction from 0{N) to 0{y/N) is due to 
the intrinsic massive parallel operation of the QC. Assuming a uniform probability distribution for the needle, for 
TV = 4 the average number of queries required by a conventional algorithm is 9/4 ||l6|,^. With Grover's QA the 
correct answer can be found in a single query p^ , p^ . 

Experimentally Grover's QA has been implemented on a 2-qubit NMR-QC for the case of a database containing 
four items ||l^,|l^. In experiments the sequences 
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Uq — X\Y \X2Y2GX1Y 1X2Y 2GX1X1Y 1X2X2Y 2^ 
Ui = XiY 1X2Y2GX1Y 1X2Y2GX iX lY 1X2X2Y2, 
U2 = X1Y1X2Y2GX1Y1X2Y2GX1X1Y1X2X2Y2, 
= XiY 1X2Y2GX lY 1X2Y2GX iXiY 1X2X2Y2, 



(24) 
(25) 
(26) 
(27) 



have been chosen to implement Grover's search algorithm. The subscript j of Uj corresponds to the position of the 
searched-for item in the database. In all four cases the input state is |00). The two-qubit operation G is defined by 



G 



( 




V 





g+i7r/4 









p+i7r/4 









g-i7r/4 



(28) 



and performs a conditional phase shift. 

On an ideal QC the QA's (p^ ) - (|2^) are by no means unique: Various alternative expressions can be written 
down by using the algebraic properties of the X's and y's. This feature has been exploited to eliminate redundant 
elementary operations |16|. On an ideal QC sequences ( |2^ ) - ( p7|) return the correct answer, i.e. the position of the 
searched-for item. This is easily verified on the QCE by selecting the elementary operations that implement an ideal 
QC. 



IV. PHYSICAL QUANTUM COMPUTER 
A. Single-qubit operations 



NMR uses SF pulses to rotate the spins. By tuning the frequency of the SF to the precession frequency of a 
particular spin in our case), the power of the applied pulse (= intensity times duration) controls how much the 
spin will rotate. The axis of the rotation is determined by the direction of the applied SF. The elementary model of 
an NMR experiment on a single spin (qubit 1 for example) subject to a constant magnetic field along the z-axis and 
a SF along the x-axis reads [|6| 



h\Sl + }v{Sl€Miujt \ |$(i)), 



(29) 



where |$) = | t)(^'| T) + I i), ~ 0)) is the initial state of the two-state system and we have set the phase 

(\>^ = 0. Substituting !$(<)) = e»*''^'5i'|*(t)) yields 



At resonance, i.e. w = /if , we find 



-h\{Sl sin cos h\t S\ smujtsmh\t) 



i^|5-(t)) = -.^(^f + 5i"sin2wt- cos2c^t) |*(t)). 

Assuming that the effects of the higher harmonic terms (i.e. the terms in sin lujt and cos liot) are small 
easily solved to give 

|*(t)) w e'*''^^i'/2|*(t = 0)), 
so that the overall action of a SF-pulse of duration r can be written as 

|$(t + r)> « e'^''i^^e*^''^^i/2|$(t)). 



(30) 

(31) 
©is 

(32) 
(33) 



Hence it follows that application of an SF-pulse of power Tft£ = tt will have the effect of rotating spin 1 by an angle 
of 7r/2 about the y-axis, as is clear by comparing (|^) with (|33|). 

In deriving (|3^), higher harmonics have been neglected, as indicated by the "~" sign. Instead of applying SF's 
along the x or y direction, one may also consider using SF's that rotate in the x-y plane. This leads to the TDSE ^] 
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hlSl + hl{Sl smut + Sf cos wt) 



(34) 



and instead of (53) we obtain 



|$(i + t)) = e*^''i^^e^^''?^?|$(0). 



(35) 



A QC contains at least two spins. If in experiments it is difficult to shield a particular spin from the SF, an 
application of an SF pulse affects not only the state of the resonant spin but changes the state of the other spins 
too (unless they are perfectly aligned along the z-axis). A general analytical, quantitative analysis of this many- 
body problem is rather difficult. We will study the limiting case in which the interaction between the spins has 
neglible impact on the time evolution of the spins during application of the SF pulse. As our numerical results (see 
below) demonstrate, this is the case that is relevant to the model system considered in the present paper and also to 
experiments 



14-17 



We consider the two-spin system described by the TDSE 



i— |$(i)) = \ hlS^ + hlSI + hi [Sf sin ujt -f cos ujt) + {S^ sin cut + cos ujt)\^ |$(t)) . 



(36) 



Substituting |$(t)) = e'^'^^^J+'^J ) |^'(t)) we obtain 



if^\^{t))=- {hl-u)Sl + {hl-u)S^i + hlSl + 'hlSl \^{t)). 



(37) 



Our aim is to rotate spin 1 about an angle ipi without affecting the state of spin 2. This can be accomplished as 
follows. First we choose 



uj = hi, 



(38) 



i.e. the frequency of the SF pulse is tuned to the resonance frequency of spin 1. Then (|37| ) can easily be integrated. 
The result is 



|$(t)) = e**''i('^i+'^2)e**''i'^ie"^^'"^^-^|$(0)). 



(39) 



where v„,™ = (0, /i^, hL_- h^). 

The third factor in (|39|) rotates spin 2 around the vector vi^2- This factor can be expressed as 



itS,r,-V„ 



1 OA ^^^^^1^ ^-ifh^^-h^^ -^h 



'n,ml 



(40) 



and we see that the SF pulse will not change the state of spin 2 if and only if the duration ti of the pulse satisfies 



ii|vi.2| = {hf - hiy + i^y 



47rni 



(41) 



where ni is a positive integer. 

The second factor in (p9|) is a special case of (Eo[). It is easy to see that if 



(42) 



the second factor in ( p9| ) will rotate spin 1 about ipi around the y-axis. Therefore, if conditions (p8|), ([4l|), and (^) are 
satisfied we can rotate spin 1 about without affecting the state of spin 2, independent of the physical realization of 
the QC. However, the first factor in (p|) can still generate a global phase shift. Although it drops out of the expression 
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of the expectation value of the qubits, in general it has to be taken into account in a QC calculation because this 
phase shift depends on the state of the spins. Adding the condition 



(43) 



where ki is a positive integer (fef > by definition), the first factor in (39) is always equal to one. Summarizing: If 
conditions (H), (|l]), (|2|), and (||) are satisfied we can rotate spin 1 about ipi without affecting the state of spin 2 
and without introducing unwanted phase shifts. 

A last constraint on the choice of the pulse parameters comes from the fact that 



h^^^h-^ , K^iK ; a^x,y,z. (44) 

Without loss of generality we will assume that < 7 < 1, in concert with the choice of parameters (||). 
Using constraint (|^ and conditions (|8|), (|l]), (|4|), and (Q we have 



and reversing the role of spin 1 and spin 2 we obtain 

(i-^)%^ + 3^(S)^-4 W 

where fci, ^2, ni, and n2 are positive integers. The angles of rotation about the y-axis can be chosen such that 
< < 27r and < ((52 < Stt. 

In general ( ^ ) or ( ^6[ ) have no solution but a good approximate solution may be obtained if 7 is a rational number 
and ki and ^2 are large. Let 7 — N/M where N and M are integers satisfying < < Af. It follows that the 



representation fci = kMN"^ and k2 — kNM^ will generate sufficiently accurate solutions of ( [45| ) and ( 46 ) if the integer 
k is chosen such that 

2fcA^M(M- A^) > I. (47) 
In terms of fc, A^, and M, the relevant physical quantities are then given by 

'-l!^ = 2kMN^ , ' (48) 

27r ' hi 2feM7V2 27r' ^ ' 

and 



27r 



2kM-^ 



"-2 
hf 



1 



^2 



2fcM3 27r ' 



(49) 



In our numerical experiments we use (Q) and (^9|) to determine the duration of the SF pulses for both the static and 
rotating SF's. In the latter case the SF pulses will be optimized in the sense that a pulse that rotates spin I (2) will 
hardly affect spin 2 (1) if fc satisfies condition (^7|). 

The assumption of a pure sinusoidal time dependence of the applied fields serves to simplify the analytical analysis 
given above. In experiment there is no good reason to stick to a simple time dependence of the pulses pq|. In general 



hlS'i + h^S^ + w{t)hl{Sl sinujt + Sf coscut) + w{t)h^{S^ sinujt + cosujt) |$(t)) 



(50) 



where w{t) can be almost any waveform. For u = hf, the formal solution of (pO|) reads 



|$(t)) = e"'*^(^?+^2')exp I I 



duw{u)h^Sf j exp^ 



du 



{h^2-K)s^+wiuMsi }\m) 



(51) 
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where exp^{. . . } denotes the time-ordered exponentiaL The introduction of a general form of 10(1) replaces condition 
©by 



exp^ <^ t j du 



and condition (H2f) becomes 



= 1, 



(52) 



tih"^ / duw{u) = ipi, (53) 
^0 

expressing the fact that the rotation angle (pi is determined by the power of the pulse only. Conditions (^) and ( |4^ ) 
remain the same. There are many forms of w{u) that will satisfy (53), so in this respect there is a lot of freedom in 
the choice of w{t). Finding the form of w{u) such that also ( |52| ) holds is a (complicated) optimization problem, in 
particular when the QC contains several qubits. 

To summarize: If conditions (|3|), (|l|), (|2|), and ^ are satisfied we can rotate spin 1 about (£1 without affec ting 
the state of spin 2 and without introducing unwanted phase shifts. In practice we may replace ( ^l[) and (^2|), by ( [s^ 
and (^) respectively. 

For our choice (|^) of the model parameters, 7=1/4 such that = 1 and M = 4. In general 7 will not be a ratio 
of two small integers but it may be approximated to any desired precision by a rational number. Let us consider 
the hypothetical case {N = 11, M = 40) such that 7 = 11/40 = 0.275. Then ^) reads 25520fc > 1 so that the 
choice k = 1 already yields an accurate solution to ( |45| ) and (^). However as tih\/2TT = 9680 and t2h\/2TT = 128000, 
rather long SF pulses are required to perform these nearly ideal, single-qubit operations. As this example shows, the 
duration of the pulses that implement accurate single-qubit operations will be determined by the representation of 7 
as a ratio of two (small) integers. 



-- 1 and M 



for accurate single-qubit operations. This implies tih\{l — 7) 



From (|8|) and (^ it follows that 
47rni or ki{M — N) — Mni. For = 1 and M = 4 we see that fci has to be a multiple of 4. This reasoning readily 
generalizes: For n spins ki{Mj — Nj) — MjUj for integer Uj and j = 2, . . . ,n. In other words, the frequencies of 
precession of each of the qubits have to be commensurate with each other. Otherwise systematic phase errors will 
be generated in the course of the computation. This conclusion does not depend on the peculiarities of the NMR 
technique: It holds in general. 



B. Two-qubit operations: CNOT gate 



As the CNOT sequence ( p^ ) has been constructed on the basis of model ( 
account for the fact that the two nuclear spins feel different static fields (see ( 



Iq), some modification is necessary to 
). In general the Hamiltonian reads 



Hnmb. = ^JSiS2 — h-iSi — h2S2- (54) 
Comparison of ( p^ with ( [5^ ) shows that the implemention of the CNOT operation requires additional rotations: 

CNOT — Y2e^''''^^^^^^^^^ e^"^'^^^^^^^ €^"^'^""■¥2, 

where we used the fact that Y2Y2 = 1. The extra phase shifts in ( |55| ) can be expressed in terms of single-qubit 
operations. The idendities 

^_,rihi-h)st ^ YiX[Yi = XiYlXi, (56) 

^^trihi^h)si = FaX^Fa, (57) 

define the single-spin rotations X'^, Y(, and X!^. 

As (^6|) and (|57) suggest, there are many different, logically equivalent sequences that implement the CNOT gate 
on an NMR QC. We have chosen to limit ourselves to the respresentations 

CNOTi = YiX[YiX!2Y2l'Y2, (58) 

CNOT2 = YiX[X'2YiY2l'Y2, (59) 

CiVOTa = XiYlX'^Y2XiI'Y2, (60) 
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where we introduced the symbol /' to represent the time evohition e,-i''^H«MR -v^jth r = — tt/J. On an ideal QC the 
sequences (|58| ) - (|6^) give identical results. On an NMR-like QC they do not because operations such as Xi and X2 
no longer commute. We will use sequences (|8|) - ( |60| ) to demonstrate the QPP. 



C. Quantum Algorithms 

On a conventional computer an algorithm is a sequence of logical operations that defines a one-to-one relation 
between the input and output data. We expect that a conventional computer always returns the correct result, 
irrespective of the input. Also a QC should have correct (input, output) relationships. In contrast to a conventional 
computer, a QC accepts as input linear superpositions of basis states and can return superpositions as well. If a 
quantum gate correctly operates on each of the basis states, it will also do so on any general linear superposition 
unless the operation generates additional phase factors that depend on the input state. Of course this does not happen 
on an ideal QC but on a realistic one it may. Above we have shown how to reduce unwanted phase errors that result 
from imperfections of the one- and two-qubit operations. 

For each realization of QC hardware, there is a one-to-one correspondence between the QA and the unitary matrix 
that transforms the superposition accordingly. A QA will operate correctly under all circumstances if the whole 
unitary matrix representing the QA is a good approximation to the ideal one. In other words, the magnitude and the 
phase of all matrix elements should be close to their ideal values. It is not sufhcient to have for example two different 
CNOT gates that operate correctly by themselves: Also the relative phases that they produce should match. For n 
qubits there are 2" (2" — 1) real numbers that specify the unitary matrix corresponding to a QA. All these numbers 
should be close to their ideal values, otherwise the QA is bound to produce wrong answers. 

Experimental realizations of QC's have not yet demonstrated that a QC can correctly compute the answer for 
inputs other than simple basis states. However, with the QC hardware currently available such a test is definitly 



within reach. The two simple QA's, (20) and ( pT| ) may be used for this purpose. 
In general on a physical QC, CNOP ^ 1 and hence {CNOTf in ^) and (|l|) 



does not reduce to one CNOT 

operation. The effect of the physical implemention of a QC on the logical operation of a QA will be most clear 
if we can distinguish errors due to faulty input data from those that are intrinsic to the physics of the qubits. 
Therefore we will provide the exact input state to the QA and compare the result returned by the QA with the exact 
answer. This procedure simplifies the analysis but does not touch the essence of the matter. We prepare I61&2) a-nd 
\singlet) — (|01) — |10))/V2 by starting from the state |00) and by performing exact rotations of the spins. 

In the case of Grover's database search algorithm, the representation of G in terms of the time evolution of ( |5^ ) 
reads 

G = e-"^i^2 = (,-irKSl^-^rh'^S^i^-^rHMMR ^ YiX'^Y 2YiX'{Y iB^ ^"^ , (61) 

where t = —i:/ J. This choice of r also fixes the angles the rotations, and through relations (^8|) and (^) also all 
parameters of the operations X'l and X2. 



V. SIMULATION 



A. Model parameters 

The parameters of model (^) for which e"*"^^ implements the EO's of the ideal QC are listed in Table |l[ On an 
NMR-like QC, the one-qubit operations can be implemented by applying SF pulses, as explained above. The two-qubit 
operation /' can be implemented by letting the system evolve in time according to Hamiltonian H^mr, given by (|5j). 
/' is the same for both an ideal or NMR-like QC. Note that the condition rJ = — tt yields T/27r — 1162790.6977, a 
fairly large number (compared to our reference hf — 1, see (|^)). 

The model parameters for the fixed and rotating SF's are determined according to the theory outlined above. 
We use the integer k to compute all free parameters and the subscript s — 2kMN'^ to label the results of the QC 
calculation. For reference we present the set of parameters corresponding to fc = 1 for QC's using fixed and rotating 



SF in Tables III and IV respectively. 
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B. Results 



As a standard test we execute all sequences on an implementation of the ideal QC (see Table ||). They all give 
the exact answers (results not shown). Furthermore the results (not shown) do not change if we put J = in all 
single-qubit operations, which is not a surprise in view of the fact that typical pulse durations are much smaller than 
1/2 1 J|. It is also necessary to rule out that the numerical results depend on the time step S used to solve the TDSE. 
The numerical error of the product formula used by QCE is proportional to It goes down by a factor of 

about one hundred if we reduce the time step by a factor of 10. Within the two-digit accuracy used to present our 
data, there is no difference between the results for 6 = 0.01 and 6 = 0.001. Hence we can be confident that we are 
solving the TDSE with a sufficiently high accuracy. 

In Table we present simulation results for QA's, QAi and QA2 defined by (^0|) and (|2^) respectively. It is 
clear that even the least accurate implementation (s = 8, fc = 1) nicely reproduces the correct answers if the input 
corresponds to one of the four basis states. The corresponding entries for QA2 seem to suggest that CNOTi is working 
well for s = 8. However the result for s = 16 (fc = 2) shows that the apparently good result for s = 8 is accidental, 
as we might have expected on the basis of criterion ( ^ ) (which in this case reads 24 ^ 1). In agreement with the 
theoretical analysis of Section IV.A the results converge to the exact ones for sufficiently large k, as shown in Table 
0. For small s, the difference in the accuracy with which QAi and QA2 give the correct answer clearly shows that in 
order for a QA to work properly, it is not sufficient to show that it correctly operates only when the input corresponds 
to one of the basis states. 

In the regime where phase errors are significant the QA's exhibit the QPP. This is exemplified in Tables ^ and 
VII where we show the results of using CNOT2 and CNOT3 instead of CNOTi. For fc < 32 there is a clear signature 
of the QPP: These logically identical QA's are sensitive to the order in which the single-qubit operations are carried 
out. 

The results presented in Tables |v|- |VIl| have been obtained using rotating SF's. As explained above, in this case 
a single-qubit operati on on qubit j exactly rotates qubit j about the specified angle (but perturbs the state of the 
other spin). In Table VIII we present simulation results obtained by using SF in the x or y direction only. Then 
the single-spin rotation on spin j no longer corresponds to the exact one. Nevertheless, as Table VIH shows, for 
sufficiently large s the results nicely converge to the correct anwers. Apparently, for a QA to compute correctly, it is 
more important to have the phase errors under control than to perform very accurate single-spin rotations. 

The very essence of QA's is the use of entangled states at some stage of the calculations. It is at this point that the 
QA is most sensitive to (accumulated) phase errors. As another illustration of this phenomenon, we present in Tables 



[X 



and g some typical results obtained by excuting Grover's database search algorithm on the same model QC's as 
those used in the examples discussed above. We find that reasonably good answers are obtained if s > 32, in concert 
with the observations based on QA's QAi and QA2. 

The results discussed above show effects of imperfections in the physical implementation of single-qubit operations. 
Thereby we assumed that J, and the static applied fields /if and /i| are fixed in time and known to very high precision. 
The Ising-model time evolution was used to perform two-qubit operations, leaving only the duration of this operation 
as a possible source for causing errors. In Table |x| we give examples of the extreme sensitivity of a QA to the precision 
with which the parameters have to be specified. Essentially we repeated the calculation of Table ^ for s = 256 but 
on purpose we made an error in the specification of the duration of /'. As Table XI shows, an error in the 8-th digit 
can have a devastating effect on the outcome of the QC calculation. This again is just another manifestation of the 
QPP but not really a surprise: During the application of /' the spins rotate around the z-axis with their resonance 
frequencies and A small deviation in t/2'k from its ideal value produces phase errors. Note however that 
the integer part of t/2t: is also essential to perform the correct conditional phase shift. Therefore, in practice it is 
necessary to specify the duration of the time evolution /' to at least 8 digits (for the case \J\/hl w 10^®). 



VI. SUMMARY 



On a physically realizable, non-ideal quantum computer, operations that manipulate one particular qubit also affect 
the state of other qubits. This may cause unwanted deviations from the ideal motion of the total system and lead to 
practical problems of programming quantum computers: An implementation of a quantum computation that works 
well on one quantum computer may fail on others. 

We have classified the various physical sources that lead to deviations. The most obvious one originates from the 
fact that other spins cannot be kept still during an operation on one particular spin. If these spins do not return to 
their original state when this operation is over, the quantum computation is unlikely to give correct answers [0| . 
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Proper optimization of the parameters that control the single-qubit operations can largely eliminate this source of 
errors. However, even if the operation gives almost exact results for all basis states, the operation is not necessarily 
perfect. That is, the operation generally yields a global phase factor which depends on the input states. Therefore, 
when such an operation is applied on a linear combination of the basis states, the relative phases of the basis states 
change, resulting in incorrect quantum computation. This is a second source for deviations from correct quantum 
computer operation. 

We have derived additional conditions on the parameters that control the singie-qubit operations and have obtained 
the conditions for reliable quantum computation. Unfortunately, these conditions cannot be satisfied simultaneously. 
However they can be satisfied to any precision by increasing the duration of the single-qubit operations. Using the 
controUed-NOT gate and Grover's search algorithm as examples, we have given concrete demonstrations of how the 
above mentioned problems arise and how they can be solved. 

At this moment, we do not know how to stabilize the quantum computation by controlling the evolution of the 
state of a closed quantum system. In a classical computer the presence of dissipation enables reliable computation. 
However, dissipation seems detrimental for quantum computer operation. Therefore, at this moment, the only option 
is to perform each operation as perfect as possible. The present paper shows how this may be done. 

The condition on the commensurability of the precession frequencies of the individual qubits leads to an increase 
of the execution time of single-qubit operations. Unless the precession frequencies of the qubits are the same to great 
precision, the execution time will grow rapidly with the number of qubits and substantially limit the speed of quantum 
computation. Therefore new techniques have to be developed to compensate for this loss in efficiency. Quantum error 
correction schemes that work well on an ideal quantum computer require many extra qubits and many additional 
operations to detect and correct errors. On a physical quantum computer however, the error-correction qubits will 
suffer from the same deficiencies as those exposed in this paper. Possibly, the clever use of dissipation processes may 
help to perform automatic error correction ]49| ]. All this puts considerable demands on the technology to fabricate 
qubits. It remains a great challenge to demonstrate that a QC of many qubits can perform a genuine computation in 
less time than a conventional computer. 
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TABLE I. Input and output states and the corresponding expectation values (a, 6) of the qubits for the CNOT operation. 



Input state 


a 


b 


Output state 


a 


b 




|00) 












00) 












|10) 




1 







11) 




1 


1 




101) 







1 




01) 







1 




lll> 




1 


1 




10) 




1 





TABLE II 


Model parameters for the elementary operations on the ideal QC. Parameters of model (3) that do not appear 


in this table 


are zero, except for the interaction J = —0.43 x 10 ^. The TDSE 


is solved 


using a time step 5/2n — 1. 








t/2tt 


hi 




hi 




hi 


hi 


hi 


Xi 




0.25 


1 



















X2 




0.25 





1 
















Fi 




0.25 








1 













Y2 




0.25 













1 








X[ 




1 


-0.4477 



















X'2 




1 





-1.4244 
















Yi 




1 








0.4477 













xi' 




1 


-0.6977 



















x^ 




1 





-1.6744 
















I 




-1/2J 
















-J/2 


-J/2 


I' 




-1/2J 
















1 


0.25 


G 




-1/2J 
















1 


0.25 



TABLE III. Model parameters of single-qubit operations on an NMR QC for the case (fc = 1, TV = 1, M = 4), see and 
(p9[). Parameters of model (3) that do not appear in this table are zero, except for the interaction J — —0.43 x 10^^ and the 
constant magnetic fields hi — 1 and /12 — 0.25. The TDSE is solved using a time step 5/2n — 0.01. 





t/2-k 


to 


hi 


hi 


hi 


hi 




8 


1.00 








-0.0625000 


-0.0156250 


X2 


128 


0.25 








-0.0156250 


-0.0039063 


Yi 


8 


1.00 


0.0625000 


0.0156250 








Y2 


128 


0.25 


0.0156250 


0.0039063 








Xi 


8 


1.00 








0.1119186 


0.0279796 


X'2 


128 


0.25 








0.0890262 


0.0222565 


Yi 


8 


1.00 


-0.1119186 


-0.0279796 








x'l 


8 


1.00 








0.1744186 


0.0436046 


x'-l 


128 


0.25 








0.1046512 


0.0261628 
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TABLE IV. Model parameters of single-qubit operations on an NMR QC using rotating SF's for the case (fc = 1, A'^ = 1, 
M = 4), see (^sh and (p9[). Parameters of model (3) that do not appear in this table are zero, except for the interaction 
J = —0.43 X 10^ and the constant magnetic fields h\ = 1 and /il = 0.25. The TDSE is solved using a time step 5/2-k — 0.01. 





t/2-k 


to 


hi 


h% 




h\ 


hi 


fv 


Ai 


8 


1.00 


-0.0312500 


-0.0078125 


-7r/2 


-0.0312500 


-0.0078125 







128 


0.25 


-0.0078125 


-0.0039063 


-7r/2 


-0.0078125 


-0.0039063 





Yi 


8 


1.00 


0.0312500 


0.0156250 





0.0312500 


0.0156250 


7r/2 


Y2 


128 


0.25 


0.0078125 


0.0039063 





0.0078125 


0.0039063 


7r/2 


X[ 


8 


1.00 


0.0559593 


0.0139898 


-7r/2 


0.0559593 


0.0139898 





X'2 


128 


0.25 


0.0445131 


0.0111283 


-7r/2 


0.0445131 


0.0111283 





Yi 


8 


1.00 


-0.0559593 


-0.0139898 





-0.0559593 


-0.0139898 


7r/2 


X'i 


8 


1.00 


0.0872093 


0.0218023 


-7r/2 


0.0872093 


0.0218023 





x'i 


128 


0.25 


0.0523256 


0.0130914 


-7r/2 


0.0523256 


0.0130914 






TABLE V. Expectation values of the two qubits as obtained by performing a sequence of five CNOT operations on a QC 
that uses rotating SF's to manipulate individual qubits. The initial states |10), |01), |11), and \singlet) = (|01) — |10))/\/2 
have been prepared by starting from the state |00) and performing exact rotations of the spins. The CNOT operations on the 
singlet state are followed by a 7r/2 rotation of spin 1 to yield a non-zero value of qubit 1. The subscripts in as and bs refer to 
the time s = r/27r = 2kMN^ and determine the duration and strength of the SF pulses through relations ( |48| ) and (|49|), see 
Tab le ffV for the example of the case s — 8. The CNOT operation itself was implemented by applying sequence CNOTi given 
by (|58[). On an ideal QC, CNOT* is the identity operation. The results obtained on an ideal QC are given by a and 6. 



Operation 


a b 


0.8 ^8 


fflie bie 


132 ^32 


ffl64 &64 


1256 ^256 


(CA'Ori)"|00) 
{CNOTiYlW) 
(CA^OTif 01) 
(CA^OTif 11) 
Yi{CNOTif\singlet) 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.90 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.03 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.58 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.88 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.99 1.00 


TABLE VL Same as Table ^ except that instead of CNOTi sequence CNOT2 given by (|59|) was used to perform the 
quantum computation. 


Operation 


b 


as. bs 


fflie bie 


132 632 


ffl64 ^64 


1256 ^256 


{CNOT2f\QQ) 
(CATOTaf 10) 
(CA^OTaf 01) 

(CAror20'^|ii) 

Yi{CNOT2f\singlet) 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 


0.24 0.76 
0.76 0.24 
0.24 0.24 
0.76 0.76 
0.98 0.24 


0.50 0.26 
0.50 0.74 
0.51 0.74 
0.50 0.26 
0.95 0.74 


0.20 0.07 
0.80 0.93 
0.20 0.93 
0.80 0.07 
0.98 0.93 


0.06 0.02 
0.95 0.98 
0.06 0.98 
0.95 0.02 
0.99 0.98 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 
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TABLE VII. Same as Table ^ except that instead of CNOTi sequence CNOT3 given by ( |60| ) was used to perform the 
quantum computation. 



Operation 


a b 


0.8 ba 




132 bz2 


ffl64 664 


0256 6256 


(CNonYioo) 

{CNOT3Y \10) 
{CNOT'iYlOl) 
(CNOTsflll) 
Yi{CNOT3f\singlet) 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 


0.23 0.76 
0.77 0.24 
0.23 0.24 
0.77 0.76 
0.79 0.24 


0.50 0.26 
0.50 0.74 
0.51 0.74 
0.50 0.26 
0.55 0.74 


0.20 0.07 
0.80 0.93 
0.20 0.93 
0.80 0.07 
0.82 0.93 


0.06 0.02 
0.95 0.98 
0.06 0.98 
0.95 0.02 
0.95 0.98 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 



TABLE VIII. Same as Ta ble M except that instead of rotating SF's, SF's along either the x or y-axis were used to manipulate 
individual qubits. See Table Ilflfor the example of the set of model parameters for s = 8. 



Operation 


a b 


08 bs 


fflie bie 


132 632 


ffl64 fe64 


0256 ^256 


(CNOT',Y\m) 
{CNOTif\\0) 
(CNOTT^f\Ql) 
[cNOT^f 
Yi{CNOTiY'\singlet) 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 


0.00 0.03 
1.00 1.00 
0.00 0.97 
1.00 0.00 
0.02 0.98 


0.00 0.01 
1.00 1.00 
0.00 0.99 
1.00 0.00 
0.45 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.17 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.70 1.00 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
0.98 1.00 



TABLE IX. Expectation values of the two qubits as obtained by running Grover's database search algoritm on a QC that 
uses rotating SF's to manipulate individual qubits. The subscripts in as and bs refer to the time s = r/27r = 2kMN^ and 
determine the duration and strength of the SF pulses through relations (^) and (^), see Table ^ for the example of the case 
s = 8. The results obtained on an ideal QC are given by a and b. 



Item position 


a 


b 


as 


bs 


016 


bie 


132 


^32 


064 


&64 


0256 


^256 





0.00 


0.00 


0.48 


0.53 


0.15 


0.16 


0.04 


0.04 


0.01 


0.01 


0.00 


0.00 


1 


1.00 


0.00 


0.52 


0.50 


0.85 


0.15 


0.96 


0.04 


0.99 


0.01 


1.00 


1.00 


2 


0.00 


1.00 


0.55 


0.48 


0.15 


0.84 


0.04 


0.96 


0.01 


0.99 


0.00 


1.00 


3 


1.00 


1.00 


0.45 


0.50 


0.85 


0.85 


0.96 


0.96 


0.99 


0.99 


1.00 


0.00 



TABLE X. Same as Table except that instead of rotating SF's, SF's along either the x or y-ax;is were used to manipulate 
individual qubits. See Table fl¥for the example of the set of model parameters for s — 8. 



Item position 


a 


b 


as 


bs 


016 


bie 


032 


^32 


064 


&64 


0256 


b256 





0.00 


0.00 


0.92 


0.91 


0.39 


0.35 


0.11 


0.10 


0.03 


0.03 


0.00 


0.00 


1 


1.00 


0.00 


0.09 


0.91 


0.61 


0.36 


0.89 


0.10 


0.97 


0.03 


1.00 


1.00 


2 


0.00 


1.00 


0.95 


0.10 


0.36 


0.65 


0.10 


0.90 


0.03 


0.98 


0.00 


1.00 


3 


1.00 


1.00 


0.05 


0.09 


0.64 


0.64 


0.90 


0.90 


0.97 


0.97 


1.00 


0.00 
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TABLE XI. Same as Table except for a change in the duration of the operation /'. (a^5'g,6^5^6): r/27r = 1162790.4977; 
(a256>&256): "^/St!- = 1162790.5977; (aatL&Se): ^/S^r = 1162790.6977 (correct value); (a^t'e'&Se): ^/S^r = 1162790.7977; 
(^sl^Se): ^/2^ = 1162790.8977. 



Operation 


a b 


"256 "256 


"256 "256 


"256 "256 


"256 "256 


"256 "256 


(Civori)='|oo) 

(CNOTiflW) 
{CNOTiflOl) 
{CNOTif\n) 
Yi{CNOTif\singlet) 


0.00 0.00 
1.00 1.00 
0.00 1.00 
1.00 0.00 
1.00 1.00 


0.00 0.52 
1.00 0.48 
0.00 0.48 
1.00 0.52 
0.99 0.50 


0.00 0.16 
1.00 0.87 
0.00 0.84 
1.00 0.13 
0.09 0.85 


0.00 0.00 
1.00 1.00 
0.00 0.00 
1.00 1.00 
0.99 1.00 


0.00 0.13 
1.00 0.84 
0.00 0.87 
1.00 0.16 
0.01 0.85 


0.00 0.48 
1.00 0.48 
0.00 0.52 
1.00 0.52 
0.99 0.50 
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